The Fundamental Physics of Shear Bands in Amorphous Solids 
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The fundamental instability responsible for the shear localization which results in shear bands in 
amorphous solids remains unknown despite enormous amount of research, both experimental and 
theoretical. As this is the main mechanism for the failure of metallic glasses, understanding the 
instability is invaluable in finding how to stabilize such materials against the tendency to shear 
localize. In this Letter we explain the mechanism for shear localization under shear, which is the 
appearance of highly correlated lines of Eshelby-like quadrupolar singularities which organize the 
non-affine plastic flow of the amorphous solid into a shear band. We prove analytically that such 
highly correlated solutions in which N quadrupoles are aligned with equal orientations are minimum 
energy states when the strain is high enough. The line lies at 45 degrees to the compressive stress. 



In Fig. [I] we show the typical failure of a sample 
of metallic glass when subjected to compressive stress. 
When the stress exceeds some critical level (known 
as the yield stress), the sample, rather than flowing 
homogeneously in a plastic flow, localizes all the shear 
in a plane that is at 45 degrees to the compressive stress 
axis, and then breaks along this plane [l[. Despite con- 
siderable amount of research [2|-|5|, the fundamental in- 
stability that gives rise to this spectacular phenomenon 
remains elusive. In this Letter we close this gap. 

Amorphous solids are obtained when a glass-former is 
cooled below the glass transition 0, @, to a state which 
on the one hand is amorphous, exhibiting liquid like or- 
ganization of the constituents (atoms, molecules or poly- 
mers), and on the other hand is a solid, reacting elasti- 
cally (revcrsibly) to small strains. There is a large vari- 
ety of experimental examples of such glassy systems, and 
theoretically there are many well studied models [8l4l0| 
based on point particles with a variety of inter-particle 
potentials that exhibit stable supercooled liquids phases 
which then solidify to an amorphous solids when cooled 
below the glass transition. Typically all these materi- 
als, both in the lab and on the computer, begin to have 
plastic (irreversible) responses when the external strains 
increases beyond some limit. All these systems also ex- 
hibit a so-called yield-stress above which the material 
fails to a plastic flow, either homogeneously or by shear 
localization as seen in Fig. [T] 
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FIG. 1: A typical example of the failure of a metallic glass 
sample when subjected to compressive stress. The material 
localizes the stress in a plane that is at 45 degrees to the 
compressive stress axis, and then breaks along this plane [lj. 



To understand this phenomenon one needs to briefly 
review recent progress in understanding plasticity in 
amorphous solids Below we deal with 2- 

dimcnsional systems composed of N point particles in an 
area A, characterized by a total energy U(ri,r2, ■ ■ ■ r n ) 
where is the position of the z'th particle. Generaliza- 
tion to 3-dimensional systems is straightforward if some- 
what technical. The fundamental plastic instability is 
most cleanly described in athermal (T = 0) and quasi- 
static (AQS) conditions when an amorphous solid is sub- 
jected to quasi-static strain, allowing to system to regain 
mechanical equilibrium after every differential stress in- 
crease. Higher temperatures and finite strain rates in- 
troduce fluctuations and lack of mechanical equilibrium 
which cloud the fundamental physics of plastic instabili- 



ties with unnecessary details [llj. We choose to develop 
the theory for the case of external shear since then the 
strain tensor is traceless, simplifying some of the the- 
oretical expressions. Applying an external shear, one 
discovers that the response of an amorphous solids to 
a small increase in the external shear strain #7 (we drop 
tensorial indices for simplicity) is composed of two con- 
tributions. The first is the affine response which simply 
follows the imposed shear, such that the particles posi- 
tions Ti = Xi,yi change via 



vi = y'i- 



(i) 



This affine response results in nonzero forces between the 
particles (in an amorphous solid) and these are relaxed 
by the non-affine response Ui which returns the system 
to mechanical equilibrium. Thus in total — > r[ + u,. 
The nonaffinc response Ui solves an exact (and model 
independent) differential equation of the form 0, llij 



where Hi, 



dui 

d'y 

d 2 U(r±,r2,--- 



(2) 



trix and H 4 = 9 U{r d ^g 2 r \ 



is the so-called Hessian ma- 
- is known as the non-affine 
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FIG. 2: (Color Online). Left panel: the localization of the 
non-affine displacement onto a quadrupolar structure which 
is modeled by an Eshelby inclusion, see right panel. Right 
panel: the displacement field associated with a single Eshelby 
circular inclusion of radius a, see text. The best fit parameters 
are a « 2.5 and e* ~ 0.1 with a Poisson ration v — 0.215. 
To remove the effect of boundary conditions, the best fit is 
generated on a smaller box of size (x, y) € [25.30, 75.92] 



force. The inverse of the Hessian matrix is evaluated af- 
ter the removal of any Goldstone modes (if they exist). A 
plastic event occurs when a nonzero eigenvalue Ap of H 
tends to zero at some strain value 7p. It was proven that 
this occurs universally via a saddle node bifurcation such 
that A p tends to zero like Ap ~ y/jp — 7 [HI . For values 
of the stress which are below the yield stress the plastic 
instability is seen Q as a localization of the eigenfunc- 
tion of H denoted as ^p which is associated with the 
eigenvalue Ap, (see Fig. [2] left panel). While at 7 = all 
the eigenfunctions associated with low-lying eigenvalues 
are delocalized, ^p localizes as 7 — > -fp (when Ap — > 0) 
on a quadrupolar structure as seen in Fig. [5for the non- 
affine displacement field, when the plastic instability is 
approached. These simple plastic instabilities involve the 
motion of a relatively small number of particles but the 
stress field that is released has a long tail. 

When the strain increases beyond some yield strain, 
the nature of the plastic instabilities can change in a fun- 
damental way [f| . The main analytic calculation that is 
reported in this Letter shows that when the stress built 
in the system is sufficiently large, instead of the eigen- 
function localizing on a single quadrupolar structure, it 
can now localize on a series of Af such structures, 
which are organized on a line that is at 45 degrees 
to the principal stress axis, with the quadrupolar 
structures having a fixed orientation relative to 
the applied shear. In Fig. [3] we show a typical stress 
vs. strain curve of a glassy material (the details of the 
simulations are presented in the accompanying material 



15j), and show the non-affine displacement field associ- 
ated with the plastic instability after exceeding the yield 
stress as indicated by the ellipse on the curve. Fig. [4] 
shows the non-affine field that is identical to the eigen- 
function which is associated with this instability, clearly 
demonstrating the series of quadrupolar structures that 
are now organizing the flow such as to localize the shear 
in a narrow strip around them. This is the fundamen- 
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FIG. 3: (Color online). The energy vs. strain in a typical nu- 
merical simulation in our system. Shown are the points on the 
curve for a regular plastic event involving a single quadrupo- 
lar structure, (marked in a red circle on the main graph and 
in insert (a) which is a blow up of the same graph), and also 
the point on the curve that results in a plastic instability lead- 
ing to a shear band, yellow ellipse on the main graph and the 
blow up in insert (b). The regular plastic event is not even 
seen without a blow up. 
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FIG. 4: (Color Online). Left panel: The nonaffine displace- 
ment field associated with a plastic instability that results in 
a shear band. Right panel: the displacement field associated 
with 7 Eshelby inclusions on a line with equal orientation. 
Note that in the left panel the quadrupoles are not precisely 
on a line as a result of the finite boundary conditions and the 
randomness. In the right panel the series of N = 7 Eshelby 
inclusions, each given by Eq. ((S]> and separated by a distance 
of 13.158, using the best fit parameters of Fig. [2] have been 
superimposed to generate the displacement field shown. 



tal shear banding instability. Note that this instability 
is reminiscent of some chainlike structure seen in liquid 
crystals, arising from the orientational elastic energy of 
the anisotropic host fluid [lij ]. and ferromagnetic chains 
of particles in strong magnetic fields 17 (. 

To explain why this mechanism for shear banding can 
appear only at values of the stress that exceed the yield 
stress, we turn now to analysis. As a first step Q 
we model the quadrupolar stress field which is associ- 
ated with the simple plastic instability as a circular Es- 
helby inclusion of radius a and a traceless eigenstrain 
e a/3 = e*(2n a np — 8 a p) where n is a unit vector |l8|. The 



inclusion is inserted in a homogeneous elastic medium 
with Young modulus E and Poisson ratio v which is 
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subjected to a homogeneous shear strain e^. It was 
shown by Eshelby that inside the inclusion we have, due 
to the effect of the constraining elastic medium, a con- 



stant strain e c „ 



c a p — S a/ 3js^* s where S a p^s 
tensor for any elliptical inclusion. For a circular inclu- 
sion the Eshelby tensor reads 



is a constant 



4i>- 1 3 -4i> 



8(1-!/) 



(3) 

Computing accordingly we find that the constrained 
strained is proportional to the eigenstrain, i.e. e c a p = 
?,7 ih \ e*n. Inside the inclusion, using the fact that a 

traceless strain field e induces a stress field <x = jq^je 
and displacement field u a = e a pXp (with X denoting an 
arbitrary cartesian point in the material) we find in the 
inclusion 



-E 



?ap ~ a/3 

Z-Av 



'a/3 



4(1 -V 2 ) 



2\ Oif) 



: *af3 + e ^]^/3 



C/3 , (4) 
(5) 



Now outside of the inclusion the displacement field can be 
written as u a (X) = u c a {x) + e^gXp where u c a (x) solves 
the bi-Laplacian equation 



V 2 V 2 <(X) = , 



(6) 



subject to continuity on the surface of the inclusion and 
zero at infinity. Remembering the radial solutions of the 
bi-Laplacian in 2-dimensions (i.e. 1, lnr, r 2 , r 2 lnr), we 
write the most general displacement field that is linear 
in the traceless eigenstrain that tends to zero at infinity: 



aj3 



Be 



Pi 



Ce 



d lnr 



<9 3 (r 2 lnr) 
^ dX a dX»dX~, 



d 3 lni 



dXadXpdX- 



(7) 



since the third derivative of r 2 vanishes identically. We 
determine the coefficients A, B, C as usual by fitting the 



boundary conditions. The calculation is lengthy but 
standard (see a line-by-line solution in the accompanying 
material) with the final result 



<(X) = 



(8) 



4(1 - v) \r 

e* /aN2 
r 



2(1-1/) 



2(1 - 2v) 
1 



a\ 2 i r2(n- X 
r J J L r 2 



2n a n ■ X — X a 
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Having the displacement field associated with each Es- 
helby inclusion at hand, we can now turn to the main 
calculation of the energy of n such inclusions arranged at 
random positions in the material and with a random ori- 
entation of their quadrupolc. Denoting the n inclusions 
with an index i = 1, 2, • ■ ■ n and the externally induced 
stress field by er^ we can write the total energy of the 
material as E = E°° + E inc + E cs h + E mat . The first is 
due to the externally induced shear interacting with the 
strain field of the inclusions, the second is the interac- 
tion of the inclusions themselves (i.e. the stress of one 
with the strain of the other) . The third is the self energy 
of the inclusions and the fourth the self energy of the 
strained material without inclusions. For the purpose of 
this calculation we need only the first two: 



E°° = - 



2 71 
i=l 



(9) 



E\ r 



na 
~2 



Here Rij is the distance between two inclusions and the 
stress field ^^(Rij) is that induced by inclusion i at a 
distance Rij away. This stress can be readily computed 
from Eq. [S] The computation of E- lnc is very lengthy, 
and is reproduced in the supplementary material in full. 
The final result is 



+4[2(1 - 2u) + (f) 2 ][2(n« • h^f - 1] - 8[1 - 2(-^) 2 ][2(n« • r t] f - l][2(n« • r l} f 1] (10) 
+32[1 - (f) 2 ][nW ■ fi 3 -nW • r tf n« ■ - (n® ■ r y ) 2 (n«» ■ r^) 2 ]} . 
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On the other hand we find Our task is now to find the configuration of M Eshelbys 

ira 2 E^e* 



E°° = -- 



4 



that minimize the total energy. Obviously, if the external Substituting this result in Eq. (flOl) simplifies it consid- 
strain 7 is sufficiently large, we need to minimize E°° crably. We find 
separately, since it is proportional to 7. The minimum 
of (fTTj) is obtained for 

n«=n« = -^. (12) 



ira 2 (e*) 2 E 
8(1 -v 2 ) 



V) tin < tMi tin 



8[l-2( — ) 2 ][2(n-r) 



+32[l-(-^) 2 ][(n.f) 2 -(n.r) 4 ]} 



I] 2 
(13) 



We can find the minimum energy very easily. Denote 
x = (n ■ r) 2 , and minimize the expression A[2x — l] 2 — 
B[x — x 2 ]. The minimum is obtained at x = 1/2, or 
cos (j) = y/l/2. We thus conclude that when the line of 
correlated quadrupole forms under shear, this line is in 
45 degrees to the compressive axis, as is indeed seen in 
experiments. Note that for other external strains which 
are not consistent with a traceless strain tensor (or in 
3-dimcnsions) this conclusion may change. 

The physical meaning of this analytic result is that 
it is cheaper (in energy) for the material to organize Af 
quadrupolar structures on a line of 45 degrees with the 
compressive stress, all having the same orientation, than 
any other arrangement of these J\f quadrupoles, includ- 
ing any random distribution. This explains why such 
a highly correlated distribution appears in the strained 
amorphous solid, and why it can only appear when the 
external strain (or the built-up stress) are high enough. 
This fact, in addition to the observation that such an ar- 
rangement of Eshelby quadrupoles organize the displace- 
ment field into a shear band, explains the origin of this 
fundamental instability. 

It should be noted that in the present calculation we 
did not predict the number of quadrupoles that appear at 
the instability. To achieve this we must consider the other 
terms in the total energy (i.e. i? es h and J3 ma t above), and 
this is beyond the scope of this Letter. This calculation 
and the resulting theoretical estimate of the yield stress 
will be presented in later publications. 
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